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Applications of the integrable system techniques to the non-equilibrium transport problems 
are discussed. We describe one-dimensional electrons tunneling through a point-like defect 
either by the s-d exchange (Kondo) mechanism, or via the resonanse level (Anderson) mecha- 
nism. These models are potential candidates to be solved exactly in the presence of arbitrary 
external bias. We draw attention also to several mesoscopical systems which can be tack- 
led by the massless form-factor approach, as perturbations of integrable models. The basic 
unperturbed model is the massless sine-Gordon model with the interaction (cosine) term 
restricted to one point, which is integrable. It is being perturbed by the second interaction 
term, which destroys integrability. Quasi-exact results can be obtained by making use of the 
basis of massless quasiparticles of the sine-Gordon model. 



INTRODUCTION 

Bethe ansatz technique (BA) is a powerfull tool for solving strongly interacting systems. 
Among its most interesting applications are the exact solutions of Kondo and Anderson 
models of impurities in the metals, sin- Gordon/ Thir ring model, etc. The central issue of 
the Bethe ansatz is a phenomenon of factorized scattering which allows for solving the the- 
ory completely starting from the two-particle scattering matrices as the only input data. 
Factorized scattering approach turned out to be more fruitful than thought earlier in the 
traditional frame of BA. Withing the standard scope of BA, as we define it, lies solving the 
bare Hamiltonian for its excitation spectrum, the distribution functions and, eventually, the 
partition function, which contains the full information on the equilibrium properties. Factor- 
ized scattering concept, however, led further to the development of new methods, providing 
us with additional important physical quantities: the series expansion of the correlation 
functions that follows from the form-factor approach and the finite-size ground state energy 
of Zamolodchikov's thermodynamic BA. The attitude which is involved here (as opposed to 
the traditional BA) is similar in spirit to the Landau theory of Fermi liquids: the funda- 
mental objects which characterize the behavior of the system are the dressed quasiparticles, 
rather than the bare particles of the free problem. Free particles, after the interactions have 
been turned on, will acquire the screening cloud of virtual pairs that will follow it, resulting 



1 The author was supported by Kreitman and by Kaufmann fellowships, and partially funded by the Israel 
Science Foundation, grant No 3/96-1. These notes were inspired by the discussions with Natan Andrei, Yuval 
Gcffcn and Yigal Meir. 



in the formation of new objects - quasiparticles that substitute the free electrons in the 
non-interacting system. The spectrum, density and interactions of quasiparticles completely 
determine the physical properties of the system. In practice it is hardly possible to deter- 
mine the properties of quasiparticles exactly, starting from the non-interacting system and 
switching on interactions, except for the special integrable cases, when the Bethe ansatz is 
applicable. Integrable quasiparticles are, however, very special as far as their interactions are 
concerned: they are stable, meaning that the conservation laws force the momentum conser- 
vation of each individual quasiparticle, and their scattering is factorized into the two-particle 
processes (in some more exotic integrable models there can in principle exist quasiparticles 
with imaginary energy - finite lifetime, so-called monstrons, but we do not consider them 
here) . Adopting integrable quasiparticle basis, one can do the best of it in describing also 
non-integrable models as perturbations of the integrable ones. This direction was considered 
in the works of Delfino, Mussardo and Simonetti [Nucl.Phys. B473 (1996) 469]. The quasi- 
particles approach turned out also to be effective in the exact description of non-equilibrium 
properties of integrable systems as shown in the works of N.Andrei [Phys.Lett. A87 (1982) 
299], Fendley, Ludwig and Saleur [Phys.Rev. B52 (1995) 8934], Peres, Sacramento and 
Carmelo ||cond-mat / 9 709 1 44]] . In the present notes we are trying to pursue further this 
approach and sketch a few possible applications of integrable quasiparticle techniques to 
the tunneling transport problems, limiting ourselves to the technically simplest cases that 
naturally lead to more realistic and complicated models. 



I. EXACT NON-LINEAR TUNNELING CURRENT 



In this section we describe two ID models of electrons tunneling through a point-like de- 
fect either by the s-d exchange (Kondo) mechanism, or via the resonanse level (Anderson) 
mechanism. These models are potential candidates to be solved exactly in the presence of 
arbitrary external bias by means of the Bethe-ansatz technique. 

1. s-d exchange model of zero-bias tunneling anomalies. 

The model we focuse on dates back to the original work of J.Appelbaum [Phys.Rev.Lett. 
17 (1966) 91]. Imagine two pieces of metal with different chemical potentials, separated 
by an insulating material. Tunneling of electrons is made possible at one special point, as 
if the two metals were connected to each other adiabatically by a "weak link." Suppose 
that, by integrating out transversal modes, one can reduce the problem to effectively one- 
dimensional, with a point-like impurity at z — 0. The precise nature of the transversal 
modes is not essencial as far as only the low-energy physics is concerned and the transitions 
between different modes can be neglected. Consider first the right hand side metal. In the 
low-energy limit one can linearize the electron spectrum near the Fermi energy, leading to 
the Hamiltonian H = J2 e ea 'Re a Re + £0>L-e a L-e- It is important that the states responsible for 
the low-energy physics can be labeled by a single parameter, e. It is convenient to introduce 

* = ty L + ty Rl z>0 (1) 
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where $ l,r are the one-dimensional left and right moving electron fields, build by the wave- 
functions of effective ID problem: ^l,r = ±J2e e±ltZettt(1 e- These fields obey the usual 
anti-commutation relations. With the boundary condition imposed, 

(0) = * L (0) + * fl (0) = 0, (2) 

situation corresponding to the full reflection. One can map one of the fields, say, right- 
moving, to the negative semi-axis preserving the boundary condition at zero as follows 

* L {-z) = -* R (z) (3) 

*l(0) = ^l(O) (4) 

The same "unfolding" procedure can be applied to the left piece of metal to map the problem 
eventually to the full axis. We end up with two species of free Dirac electrons on the line, 
to which we assign the flavor index m = 1,2. The Hamiltonian describing this problem is: 

/oo /*oo 
y+ m d z V am dz + Urn f g + ra f„ (5) 

-OO in J — OO 

a,m=i,z <r,rn=l,2 

where we introduced the bias term which accounts for the fact that the left and right metals 
are at different chemical potentials /ii i2 . We set \i\ = — /x 2 = £t an d rewrite the last term as 

/oo 
- tf^flV) (6) 

The Hamiltonian (§) possesses SU(2) flavor and SU(2) spm symmetries, as well as U(l) charge 
symmetry, and the total number of type 1 and type 2 particles is conserved separately. Now 
we introduce a weak tunneling between the metals via an impurity spin S — 1/2 at the point 
z = 0: 

Hint = EE J?*ta(0)<&>*M0)S x (7) 

A=l ab 

where a x are Pauli matrices and Jf' are coupling constants. We make a natural assumption 
that Jf = Jf and J 12 = J 21 , J 11 = J 22 . The full Hamiltonian H +H int resembles the two- 
channel Kondo model, with however additional intra-channel interaction. The interaction 
term J 12 breaks the SU(2)f lavor symmetry down to its JJ(l)^ lavor subgroup. As a resut, the 
number of type 1 and type 2 particles is no longer conserved separately; however, the total 
number of type 1 and type 2 particles is still a good quantum number: 

gcharge = J ^ + ^ + (g) 

The Sl lavor charge corresponding to U(l)^ lavor subgroup also remains conserved: 

s flavor W +7^*6 = / #1~*2 + (9) 
ab 

Here t x is the first Pauli matrix acting in the flavor space, and the summation over the spin 
index is implicitly assumed. 
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We are not going to discuss the experimental relevance of this particular model, i.e. 
what physical systems it actually describes. We believe that it is one of the simplest and 
yet non-trivial theoretical models, and therefore it is instructive to learn as much as possible 
from it. We shall mention, however, that more complicated and realistic models, such as 
the one introduced by Vladar and Zawadowski [Phys.Rev. B28 (1983) 1564], are its direct 
generalizations as far as the methodology is concerned. The reader interested in the relevant 
experimental advances can see the papers of Ralph and Buhrman [Phys.Rev. Lett. 69 (1992) 
2118, Phys.Rev. B51 (1995) 3554]. Some exact theoretical results in this field, based on 
the conformal field theory and the Hershfield's formalism, include the exact scaling func- 
tions found by Delft, Ludwig and Ambegaokar [ |cond-mat/970204TS| , submitted to Annals of 



Physics]. It is worth mentioning that the one-dimensional approximation that we adopted 
may not describe correctly certain microscopic quantities, such as the local density of states 
in the realistic point contacts, as was noticed by Ulreich and Zwerger ||cond-mat/9710i74] , 



submitted to Europhys.Lett.]. However, it does work extremely well for the linear conduc- 
tance at zero temperature [Phys.Rev. B50 (1994) 17320]. 

One can bosonize the fermions according to the standard procedure and express the 
Hamiltonian in terms of four left-moving boson fields: 

tf+ m =: e**™ : (10) 

Following Emery and Kivelson [Phys.Rev.B47 (1992) 10812], introduce linear combinations 
of the bosons as follows: ^ 

<t>c= 2 + ^1 + ^2 + ^2) (11) 
<Ps = 2(^1-^1 + ^2-^2) (12) 
0/ = -(<p n + <pu -^ T 2 -fist) (13) 

4>sf = - fa - vw + vi?) (14) 

The field C represents the total charge degree of freedom in two channels and decouples 
from the Hamiltonian. The resulting Hamiltonian depends on the remaining three fields as 
follows: 

#o = |J E J(d z <p m f + ^J d-A f (15) 

m=s,f,sf '' 

H int = A[Jl 2 cos 0/(0) + Jl 1 cos0 s/ (O)] • [cos 0,(0) S x + sin0 s (O)S y ] (16) 

7 11 

- 4J 2 12 sin 0/(0) sin0 s/ (O)S 2 + -*-d g <l> a {0)S z (17) 

7T 

Interaction term ([16]) resembles the one-channel Kondo model where the coupling constant 
J x became a dynamical degree of freedom, J x — > [J* 2 cos 0/(0) + J^ 1 cos0 s /(O)]. There exists 
a special point in the space of couplings (J^ 1 — and J^ 2 = 0), an analogue of the 

Toulouse point in the Kondo model, where the full Hamiltonian H = H Q + H int becomes 
quadratic in terms of new fermion operators and can be solved exactly. This was noticed by 
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Schiller and Hershfield [Phys.Rev.B51 (1995) 12896], who found the non-equilibrium current 
at this "Toulouse point." 

In the absense of bias, the full interacting Hamiltonian is integrable. To see this, 

let us make a canonical transformation to another basis of fermions: 

a«x = -^(*i CT + * 2ff ) (18) 

In this new basis the matrix of couplings J ab becomes diagonal with the eigenvalues = 
J 11 ± J 12 . The resulting Hamiltonian is the channel-anisotropic two-channel Kondo model, 
which was solved by means of Bethe ansatz by N.Andrei and A.Jerez [Phys.Rev.Lett.74 
(1995) 4507]: 

H = tv F J2[ (aZd z a a + btd z b a ) + £[# a+(0)<7£^(0) + J;bt(0)a x aa ,b a ,(0)]S x (20) 

& A 

y = /iE/ a X + ^X (2i) 

The advantage of this new basis is that the number of a-particles as well as b-particles is a 
conserved quantity, and the Bethe wave functions can be constructed in the sector of fixed 
quantum numbers N a and Nj,. These conserved charges correspond to the total particle 

number Q char 9e and gflavor in the old bagig of ^ ({/(^charge and JJ ^flavor symme t r i eS ). 

Operator Y , being a generator of the SU(2)* lavor symmetry, does not commute with the 
Hamiltonian, except for the channel-isotropic point J + = J~ , corresponding to J 12 = 0, 
where the SU{2)^ lavor symmetry is restored (note that in general J + > J - ). The non- 
conservation of SU(2)-f lavor charge, in particular of its component Yq, is responsible for the 
transport properties we are interested in. 

Following Andrei and Jerez, let us study one and two-particle sectors of the first-quantized 
analogue of Hamiltonian ( p0|) . Naively, the scattering matrix on impurity would be 

R a a = exp(iJ + aS), R b b = exp(U~&S), R b a = 

where the scattering is obviously diagonal in the flavor space. Respectivly, the scattering 
matrix on impurity dictates us what the scattering matrix in the bulk would be, since the 
basis of the wave-functions in the bulk is not unique and must be chosen in the consistent 
way with the impurity. Let us clarify this point. Consider one of the possible representations 
of the two-particle wave function for free ID fermions with the linear spectrum: 

iP(x 1 ,x 2 ) = 9( Xl - x 2 )(A aia2 e iklXl+ik ^ + B aia2 e iklX2+ik2X ') 

+6(x 2 - x 1 )((SA) tTlt72 e ik ^ +ik23:2 + {SB) aia2 e iklX2+ik2X ') 

Here A and B are spinors, S is for the time being arbitrary matrix acting on the spinors, 
and 6{x) is a step function. It is easy to check that ip( x i, x 2) is a solution of Schrodinger 
equation 

(-id Xl - id X2 )ip = Eip 
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with E — ki + k2- In addition, one must require tp(x\,X2) to be antisymmetric with respect 
to the simultaneous permutation of spin and coordinate indexes. This leads to the relation 
between A and B, 

B = -PSA 

and the constraint for S: 

S12S21 = 1 (22) 

where P is the permutation operator acting on spinors and S21 = PSvzP- The condition 
(P^) leaves some freedom for us to choose the matrix S (for example, S = I and S = P are 
valid solutions to (0)). The consistency between impirity and bulk scattering is represented 
by the boundary Yang-Baxter equation, 

S12R2R1 = R1R2S12 

which imposes a further constraint on S, or in other words selects bulk scattering in a 
consistent way with the impurity. With the impurity scattering matrix given above, the 
Yang-Baxter equation would be solved by the bulk S-matrix 

g pspin pflavor 

where P are exchange operators acting in spin or flavor spaces of two particles. As noticed 
by Wiegmann and Tsvelik, as well as by Andrei and Destri, such a choice of scattering 
matrices leads to wrong results. The Hamiltonian must be regularized appropriately from 
the very beginning, which modifies S-matrices in such a way that they become momentum- 
dependent (it was not possible before regularization because J is a dimensionless coupling 
and there is no scale in the problem; after the regularization the model acquires a scale 
which enters S-matrices in the dimensionless combination with momentum). Andrei and 
Destri [Phys.Rev.Lett 52 (1984) 364] suggested an elegant regularization compatible with 
the integrability: one introduces the second-order derivative term coupled to a mass scale A 
in the first-quantized Hamiltonian, A -1 <9 2 , thus "delinearizing" the spectrum. After such a 
regularization the bare scattering matrices become: 

^ /A ;;;:^ 1/2) 

(impurity scattering is still flavor-diagonal) 

q s P in (h , v _ (fcl/A + l)/J± - (k 2 /A + l)/J± - iP*" f 

bl2 {ku k2) ~ (h/A + 1)/ J± - (k 2 /A + 1)/ J± - i {2b) 
where J + (J _ ) is chosen depending on the flavor of the first and second particles, 

Sfr-ih-h) = l sinh[K{ki Sm ; 2)/A - -j [r, ® r x + r y ® r y ] + 1 [1 + r z ® r z ] 

1 sinh[/t(fci - k 2 )/A] 
+ 2sinhK fcl -A; 2 )/A + H [ ^ A ( } 
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Here r are the Pauli matrices acting on the flavor space, k and v are functions of J^. In the 
scaling limit we have v/k = J~ . The complete bulk scattering matrix is a tensor product of 
spin and flavor terms, S12 = Sy^ n ® S{^ vm . The S^ n term reflects the SU{2) spin invariance, 
while the flavor component S^ 07 " reflects the breaking of the SU(2)f lavor symmetry to 
U{l)f lavor subgroup in the presence of anisotropy. 

What is the form of the above S-matrices in the original basis of \I/-fields? Since our 
transformation fll8|)-(|l9D mixes up the flavors (and not the spins), we will list here only two 
of the above terms mapped to the original basis: 

R\ = Rl = \{R a a + R b b ) (27) 

R\ = R\ = ^(K ~ R b) ( 28 ) 

The impurity scattering is obviously non-diagonal, allowing the processes of flavor violation 
1 — ► 2, which accounts for the non-trivial transport through impurity from the left to the 
right metal. Further, the bulk scattering becomes 

§( > VOr = \ sinhKA, -T 2 )/A + iu] [Tz ® Tz + Ty ® Ty] + l [1 + Tx ^ Tx] 

1 sinh[K(fci-fc 2 )/A] r , , , 

+ ^sinhK^-^/A + H' X ® J (M) 

The non-conservation of flavor is also explicitly present in this term: the in-state |1) ® |1) 
will be mapped to the supperposition of |1) ® |1) and |2) ® |2). So, the scattering of particles 
of first kind into particles of the second kind occurs not only on the impurity, but also in 
the bulk. 



2. Non- equilibrium transport through a resonance level. 

Anderson proposed a more generic model than Kondo-type model described in the previous 
section. All the arguments of the previous section up to formula (0) should be recalled as 
they are unchanged here. Following Anderson, one assumes the existence of localized states 
(either impurity or surface states) on the interface between two metals. Such states, or for 
our purposes just one state is introduced formally as d+\0), and the coupling to bulk electrons 
takes a form (instead of equation (0)) 

H int = VE[*™(°K + #*«m(0)] (30) 

(7,771 

Besides, one adds to Hamiltonian the part which describes the localized level: 

H d = e d Y, d + a d a + dtd a d\d- a (31) 

The stationary one-particle eigenstate of Hq + Hi nt + is a superposition of three terms, 
M = J dx[ gi (x)^t(x)+g 2 (x)^t(x)}\0)+edM (32) 
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(the particle can be either on the localized level or in the bulk in one of two flavors with spin 
preserved). From it we infer the impurity scattering matrix 

R\ = R\ = R l = R l = jHE) (33) 
1 V 2 

A(E) = -2 tan" 1 (34) 

E-e d 

The Anderson model has U(l)^ lavor symmetry, like the Kondo model of the previous 
section. Notice, however, that the flavor SU (2)-symmetric two-channel Anderson model can 
be used to get the above model fl30|)-(|3lD. The S'[/(2)-symmetric Anderson model differs by 
the presence of two localized levels, di a \0) and d^O), each one coupled to its channel, 

H int = VYX^tmW.m + dL*am(0)] (35) 
<r,m 

Hd = £d d^dma- — — d^d^^dmio-dmo-i (36) 

m,d mm! a a 1 

In the 5'L r (2)-symmetric model the number of particles of each kind is conserved separately 
(particles in the bulk and on the localized level must be added together to obtain a conserved 
charge), and moreover, the SU(2) symmetry in the flavor space is present. To break this 
symmetry down to U(l) and to destroy the conservation of each kind of particles separately, 
thus reproducing the original Anderson model (j30j)-(|31~|), we introduce an additional term, 
describing hopping between localized levels: 

#12 = -tE«^ + dtJia) (37) 

a 

Then one must take the limit t — > oo to reproduce the model of interest (intuitevely, this 
limit means identifying di with d<i and calling both of them by d)^\ 

The same discussion about conserved charges as in the previous section applies to An- 
derson model. In particular, after the same canonical transformation of fields (|I~8|)-([T9l) we 
obtain the Anderson model with one kind of particles coupled to impurity, while another 
kind of particles being free. These two kinds of particles are, however, coupled to each other 
by the bias term: 

H a = iv F J2 / atd^ + VYsKWda + d+a^+e^dtda + UYsdtdvdtj-a (38) 

a a a a 

H b = tv F J2 I Kd z K (39) 

y = /iE / 4K + Ka a (40) 



2 When hopping amplitude t is very large, it is worth to work in the basis of mixed states e\ = d\ — 
c?2, e-2 = di + ^2, where the hopping and €d -terms become diagonal: (ed + £)e^ei + (e<j — £)e 2 f e2- In the 
limit t — ► oo one of the mixed localized levels (ei) becomes so high in energy that the transitions to this 
level are unlikely to take place and this level decouples. Then, after rescaling ed one can identify e2 = d (the 
spin index is assumed everywhere). 
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3. Open questions. 



Computing the non-equilibrium current implies an evaluation of the average 

d , d 



e— 
dt 



(Y ) = en-Y, f(4K + bta a ) (41) 



on the Bethe states. Fendley, Ludwig and Saleur [Phys.Rev. B52 (1995) 8934] suggested that 
at least for certain integrable models (in their case quantum Hall bar with a constriction) 
the following procedure can be used. First, based on the Bethe ansatz solution in the 
absense of bias, one identifies bulk excitations (physical quasiparticles) that carry a non- 
trivial charge/flavor. Further, using standard techniques, one calculates scattering matrices 
in the bulk and on the impurity for the physical quasiparticles and, in the absense of impurity 
but in the presence of bias, densities of states. Then, based on the probabilistic arguments, 
one employes the Bolzmann rate equation to obtain the current: 

deMe,/i)-n 2 (e,/i)]|£ 12 (e)| 2 (42) 

Note that this procedure has not been shown rigorously to coincide with (|4TD, and remains 
as a conjecture [ |cond-mat / 9708 TEQ 9710205]. Working in the original \l/-basis, we face the 
problem of non-diagonal bulk scattering - the scattering with the flavor violation 1 — ► 2 
happens everywhere, and not only on the impurity, as it was in the case of the Quantum 
Hall problem with a constriction. The flavor z-projection Yq is not diagonal on the Bethe 
ansatz states. This leads to some technical difficulties, as well as to the question of the 
applicability of ( f42|) for such situations. Equivalently, the difficulty can be understood in 
the representation (p0|) of a and b-particles and Eq (f41~|) . Operator Y (x-component of the 
flavor in this representation) strongly disturbs the Bethe state and does not allow us to work 
only with a few significant low-lying excitations. The state obtained after applying Y to the 
Bethe vacuum state involves infinite amount of excitations in the thermodynamic limit. The 
situation is rather similar to the XXZ chain and the S x operator, which perturbs spins on all 
sites simultaneously, while the Bethe states are constructed as local spin-flips with respect 
to some reference state. Thus, one says that S x excites an infinite amount of elementary 
particles. 

To summurize, the questions to be answered are: what are the physical excitations in 
the Bethe ansatz solution of ©-(0) and their scattering matrices? What are the matrix 
elements of the fields between the excited states (form-factors for the two-channel Kondo 
model)? Can one apply Eq (|42|) to the model (H)-(0)? 



II. COULOMB BLOCKADE AND THE INTERFERENCE 



We describe several problems of interest to condensed matter physics which could be solved 
by the massless form-factors approach, as perturbations of integrable models. The basic 
unperturbed model is the massless sine-Gordon model with the interaction (cosine) term 
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restricted to one point, which is integrable. It is being perturbed by the second interaction 
term, V(0), which destroys the integrability: 

1 roo 

H=- [7r 2 (x) + (a a; 0) 2 ]dx + /icos/30(O) + AK(0) (43) 

2 jo 

Quasi-exact results can be obtained by making use of the basis of massless quasiparticles of 
the sine-Gordon model. Namely, one expands the evolution operator, S = Te lX -f in the 
various matrix elements of interest, 

(TO x 2 ...) = (TOiOz . . .) + %\{(IO x 2 ...f V(<t>)) + ... (44) 

to generate a perturbation series for the physical quantities of interest. Further, for the eval- 
uation of averages in ([44]) one inserts the full set of intermediate states - massless integrable 
quasiparticles of the unperturbed sine-Gordon model: 

(O 1 O 2 ...yV(0)>= 2 (0\O 1 \n)(n\O 2 \m)...(l\fv(<f>)\0) (45) 

n,m,l... 

The form-factors (n\0\m) are known exactly. While generating the perturbation series, the 
following normalization conditions should be respected: 

<0|S|0> = (0|0) (46) 
= (1|1'> (47) 

The first equation fllB]) takes into account the renormalization of the ground state energy 
(disconnected diagrams), while the second one fl4~7|) controls the norm of one-particle states. 



1. Quantum dot 

Matveev [Phys. Rev. B51 (1995) p. 1743] argued that the physics of Coulomb blockade in the 
quantum dots can be adequately described by the following Hamiltonian, H = H + Hc + H': 



H = - 

2 



dx 



7T 2 (x) + (d x U) 



H' - 



--EcH0)-N} 2 
-Vcos[27ra(0)] 



(48) 

(49) 
(50) 



where Eq = t^t is the charging energy characterizing the dot, parameter iV is proportional 
to the gate voltage which is adjustable. Electron tunneling into the dot leads to the increase 
of energy by the amount Eq. Therefore, the tunneling conductance is suppresed. The gate 
voltage allows to control electrostatic energy of the dot: 



Er 



(Q - eNf 

2Cn 



(51) 
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where Q is the total chatge of the dot. At the values of the gate voltage corresponding 
toiV = ri + l/2 the energies of states with charges en and ein + 1) are equal, leading to 
the favorable conditions in the electron transport through the dot. Therefore, one observes 
periodic resonanses of the conducatance as a function of the gate voltage. The term H' de- 
scribes quantum tunneling (backscattering) effects, while the term H c describes electrostatic 
charging effects. The Coulomb blockade shows up in the oscillations of the ground state en- 
ergy or the average charge (Q) as a function of N, which can be calculated. Experimentally 
the capacitance C = d 2 E/dN 2 can be measured. At the perfect transmission, V = 0, the 
ground state energy is iV-independent, as can be seen from (pjf)-fl4^) by the change of vari- 
ables u — > u + N. The Coulomb blockade is completely suppressed. On the contrary, for the 
perfect reflection, V — ► <x>, the dependence on N is the most drastic, for the field is pinned 
at zero by the value n corresponding to one of the equivalent minima, and a departure of N 
from integer n leads to the increase of energy given by the term He, thus lifting the degener- 
acy. At the values N = n + 1/2 two neighboring minima still remain degenerate, leading to 
the drastic changes in the transport properties. Matveev in his work calculates perturbative 
corrections to the ground state energy. He treats the term H' as a small perturbation and 
uses the exact eigenstates of Hq + He- The form-factors of sine-Gordon model allow, on 
the contrary, treat exactly the H' term, starting from the integrable model Hq + H' and 
to perturb it by He term. The required form-factors of the field u, originating from the 
expression for the ground state energy correction, (0|Texp(i / dtH c )\0), are available. 



2. Quantum Hall interferometer. 



Chamon et al | |cond- mat/9607195] ] discussed the device based on the quantum Hall bar with 
two constrictions which allows observation of interesting interference effects. The underlying 
Hamiltonian is H = H + H tun : 



Hn 



1 

8^ 



dx 



n 2 (x) + (d x <f>y 



(52) 



Htun = 7i cos#(0) + 7 2 cos[6(d) + $ ] (53) 
where the two fields 9, <j) are related to the left-right moving fields 4>l, <pR as follows: 



V2<j) 
V28 



(54) 
(55) 



The quantity of interest is the tunelling current It($o, d) as a function of flux $ and distance 
d. For the current It it is sufficient to take its linear response limit, given by Kubo's formula. 
The current-current correlators can be obtained using form-factors method and treating one 
of the cosine terms in (|53"D as a perturbation of integrable sine-Gordon model. 



3. Double barrier in the Luttinger liquid. 
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Kane and Fisher [Phys.Rev. B46 (1992) 15233] studied various tunneling effects through the 
barriers in Luttinger liquid. The double barrier action reads: 

S = S + V J dT{cos[2^e(T,x = 0)] + cos[2^6(T,x = d) + k F d}} 

+ -$L [drW(T, x = d)-6(r, x = 0)1 (56) 

\/7T J 



where Sq is the action of the pure Luttinger liquid. In the last term the total number of 
particles between two barriers is coupled to the chemical potential of the island, Vq- In 
analogy with the problem of quantum dot above, one can adjust Vq to have the energy cost 
of adding another particle vanish. Then, one expects a resonant transmission. The mass 
scale M = k F j gd measures the fluctuations of the charge on the island: 

k F d V G 

(n) = 57 

2n M K J 

The first term in the above formula is just the background density. The resonanse is achieved 
when (n) =half- integer. In the strong barrier limit, V ^> M, the total charge of the island 
prefers to be a particular integer times e, since the discreteness of the electron charge is 
important. The Coulomb blockade supresses the transport through the island. When (n) is 
tuned to be a half-integer, an additional symmetry is present which leads to the degenerate 
states of different charge. In the opposite, weak barrier limit, V <C M, following Kane and 
Fisher one can integrate out small fluctuations of the charge around the value (n) and obtain 
an effective action: 

r bV 2 r 

S = S + aV J dTcos[2V7T0(O)] + —J dr cos [V5r0(O)] (58) 

where a, b are dimensionless constants. The second cosine term in the last formula is ir- 
relevant for g > 1/4, while the first one is relevant. At small temperatures the effective 
barrier strength a grows, making the conductance to vanish at T = 0. As for the effective 
coupling b, it decreases, leading thus to the perfect conductance. Thus, by fine-tuning the 
single parameter a = (corresponding to the vanishing of the 2kp Fourier component of 
the double-barrier scattering potential V(x)), one achieves the resonance. Qualitatively, the 
physical meaning of the resonance is the following: forbidding the one-electron backscatter- 
ing processes (V(2k F ) = 0), one creates favorable conditions for the coherent transport of 
electron pairs (V"(4fcp) ^0). It is interesting to find the conductance as a function of a, b by 
the form-factor method, treating one of the cosine terms in (|58|) as perturbation. 



4. Potential difficulties. 

In all the three problems one has to deal with IR divergences in the form-factor expansion, 
since the perturbing terms are either cos or field itself. Such divergences are peculiar for 
the masssless theories. Lesage and Saleur | |cond-mat / 9 60 81 12] ] have discussed already the 
situation. The partial resummation of the most divergent terms of the perturbation series is 
necessary. 
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In the first problem it is not clear that the perturbation theory in Ec will describe 
correctly the interesting regime of large Ec, although the rapid convergence of the form- 
factors series leaves us some hope that one is allowed to reach large enough values of E c 
while the terms of the order of E c and higher still can be neglected. 

In the second problem the usual trick of changing the basis of fields to the odd and even 
fields does not work with two constrictions. The odd and even fields are not decoupled any 
longer. 

In the third problem one has to deal with the irrelevant perturbation, for which the form- 
factors are not quite known. One can in principle overcome this difficulty by an analytic 
continuation from the space of relevant perturbations. The breather form-factors can be 
analytically continued from the sinh-Gordon model, while kink and anti-kink form factors 
can be obtained by continuing analytically the recent results of Lukyanov | |hep-th/9703190 
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